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ABSTRACT: We consider the probability P(i) that a given site remains unvisited by 
any of a set of random walkers in d dimensions undergoing the reaction A + A — > when 
they meet. We find that asymptotically P(t) ~ with a universal exponent 9 = ^ — 0(e) 
for d = 2 — e, while, for d > 2, 9 is no n- universal and depends on the reaction rate. The 
analysis, which uses field-theoretic renormalisation group methods, is also applied to the 
reaction k A — > with k > 2. In this case, a stretched exponential behaviour is found for 
all d > 1, except in the case fc = 3, d = 1, where P(t) ~ e -const.(lni) 3 / 2 _ 



In a recent letter, Derrida, Bray and Godreche 1 have found new non-trivial and ap- 
parently universal exponents associated with the zero-temperature relaxational dynamics 
of the one-dimensional Ising and Potts models. The occurrence of such exponents is sur- 
prising given the trivial nature of the conventional static and dynamic exponents in these 
models. However, the quantity considered by these authors, namely the probability P(t) 
that, starting from a random initial configuration, a given site has not been crossed by a 
domain wall, is not simply related to the usual response functions, and might be expected 
to show more interesting behaviour. The fact that a simple universal power law P(t) ~ 
is obtained for large t, however, requires some explanation. 

In this letter, we provide such an explanation within the context of a generalisation of 
this problem to arbitrary dimensionality d. Since the motion of domain walls for d > 1 is 
very difficult to treat analytically, instead we observe that, in one dimension, the motion 
and annihilation of Ising domain walls at zero temperature is equivalent to a reaction- 
diffusion process of point particles A undergoing the irreversible reaction A + A — > 0. The 
study of this problem is readily capable of generalisation to arbitrary d, and many of its 
features have already been elucidated using a number of approaches. 2 In particular, it is 
found that there is an upper critical dimension d c = 2 above which the mean density n(t) 
behaves as l/(Xt), where A is the reaction rate, as predicted by a simple rate equation 
neglecting correlation effects, while for d < 2 these effects cannot be ignored and the 
behaviour is modified to t~^/ 2 , with an amplitude independent of A. Recently, a systematic 
field-theoretic renormalisation group approach to this problem has been developed, 3 which 
not only yields the exponents but also correlation functions and universal amplitudes 
within a e-expansion. It is also straightforward to generalise the analysis to the reaction 
kA — > 0. In this case, above the upper critical dimension d c (k) = 2/(k — 1) one finds 
n(t) ~ l/(A£) rfc ( fc )/ 2 , while for d < d c , n(t) ~ t~ d / 2 with a universal amplitude. 

Within this type of reaction-diffusion problem, then, we ask the following question: 
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from a random initial condition (mean density no) at time t = 0, what is the late time 
dependence of the probability P(t) that a given site has never been visited by a walker? 
A simple approach to this problem is to note that P(t) — P(t + 5t) is the probability of a 
finding a walker at the given site (the origin, say), in the time interval (t, t + St), given that 
the origin has never been visited in the past. This will happen only if a particle happens 
to lie close to the origin at time t. Thus 

-P'{t)8t ~ DStP(t)h(t) 

where D is the diffusion constant, and n(t) is the density at a site adjacent to the origin, 
given that the origin is never visited, that is, with a repulsive potential there. 

For d > 2, only a finite fraction of particles near the origin has ever visited the origin 
in the past, so that n(t) oc n(t) ~ t- rf c(fc)/2_ For k = 2j this leads to P(t) ~ t- const -/ A , 
that is, a power law with a non-universal exponent, while for k > 2 we obtain a stretched 

4-(k— 2) /(&— 1) 

exponential behaviour P(t) ~ e -const.t p or ^ < 2, however, almost all particles 

near the origin have visited it at some time in the past, so that h(t) <C n(t). When 
2 > d > d c (k) (which is possible when k > 2), particle correlations may be neglected and 
we may use the inhomogeneous rate equation 

dn/dt = DV 2 n - k\n k , (1) 

selecting the required events by imposing the condition n(r = 0, t) = 0. This problem has 
the radially symmetric scaling solution n(r,t) = t^ dc ^l 2 f(r /(Dt) 1 / 2 ). For r 

(^) 1/2 , 

the nonlinear term is unimportant (corresponding to the fact that the density is so low 
that annihilation events rarely occur), and / satisfies Laplace's equation, with solution 
/ ~ (r/(Dt) 1 / 2 ) e , where e = 2 - d. Thus n(t) ~ n(i)t~ e / 2 , and the stretched exponential 
becomes P(t) ~ e - consti(d " dc(fe))/2 . For d < d c (k), if we assume that n(t) is suppressed 
relative to the bulk density by this same factor, we find that n(t) ~ £-<V2-e/2 = t -l ^ 
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resulting in a power law P(t) ~ ', consistent with the result of Derrida et al.. 1 However, 
to justify this argument and to demonstrate the universality of 6>, it is necessary to proceed 
more systematically. 

We first relate P(t) to an appropriate correlation function in the field-theoretic descrip- 
tion of the problem. We follow the notation and formalism of Ref. 3. Following Doi 4 and 
Peliti 5 , the master equation for the reaction-diffusion problem is encoded in a hamiltonian, 
or liouvillean, which may expressed in the 'second-quantised' form 

H = (D/b 2 ) £(aj - a])(a t - a-) - a£(1 - {a\) k )a k 

n.n. i 

where the first term is a sum over nearest neighbours, and represents a continuous-time 
random walk on a lattice with spacing 6, and the second term represents the annihilation 
process kA — > 0. The time translation operator is e~ Ht , which may be written as a path- 
integral by dividing the interval (0, t) into small slices of duration At. At each slice a 
complete set of coherent states 



J e -#^e&°V><0|e# a #;F#i 



is inserted (lattice labels are suppressed for clarity). The matrix elements 
(Ole^+AtV^lO) = e#+At& then give rise, when combined with the measure factors 
e -<t>t<f>t^ t Q the time-derivative piece in the action 

S = J (<p(t)*d t <p(t) + L»V0* • V0 + \(<p* k - l)(p k )dt 

in the limit At — > 0. 

In the second quantised formalism, the probability P(t) that a given site, which may 

chosen to be the origin j = 0, is never visited is simply given by inserting the operator 

Sx at each time slice. Thus, at the origin, instead of the matrix element given above, 
a' a ,0 

one should compute 

(Ole^+At^, e <M| ) = i 

X 1 aU,0 1 1 

4 



which means that, to leading order in At, the factor from the measure is not cancelled, 
corresponding to an insertion of [j t e~^ ' f ^ 0,f in the path integral. In the continuum limit 
the lattice fields are rescaled so that 0*0 — > 6 rf 0*0, so that finally 

P ( t ) = J V(j)*V(j)e~ h I ^(0,t)<f,(0,t)dt e -S[<p*,<P) 

where h = b d /At. 

As explained in Refs. 3,5, in order to compute statistical averages with respect to 
the action S, it is necessary to evaluate projections onto the state (0| Ylj eaj - ^ ^ s then 
convenient to make the shift 0* = 1 + 0, so that annihilates this state when acting to 
the left. The required insertion then has two pieces 

e -h J* 0(O/)0(O,f )df e -h J* <f>(0,t')dt> 

In the first factor, the time integral can be taken up to infinity, and thus the term 
h J 0(0, t')(f)(0, t')dt' may be regarded as a repulsive potential and included in the action 5, 
while the second factor is the piece whose expectation value we now wish to evaluate with 
respect to this modified action. It is convenient to rewrite this as a cumulant expansion 

P(t) ~ ( e - h Io^ ^ dt ) = e -^/o^(0 ) t')>* / +^/o7o'^( '*')^(0/')>c^^"+- (2) 

Dimensional analysis then dictates that h has dimension (wave number) 2 "^, so the 
additional interaction term is irrelevant for d > 2, and may therefore be neglected in 
studying the late time asymptotics. For d > d c (k), the arguments of Ref. 3 show that all 
loop corrections to the field theory are also irrelevant. The sum of the tree diagrams is 
then given by the solution to the naive rate equation, (0(0, t)) ~ l/(\t) dc ( k ^ 2 . (In fact 
the amplitude will be modified in the neighbourhood of the origin, but not the exponent, 
for d > 2.) On substitution into the first term of the cumulant expansion, this then leads 
to the same result as the earlier naive argument. The higher order terms in the cumulant 

5 



expansion all involve at least one more power of A, and hence their integrals are down by 
successive powers of £-( rf -4(fc))/2_ The appearance of the borderline dimensionality d = 2 
is simply related to the recurrence property of random walks. 

For d < 2, however, the h interaction is either marginal or relevant, and it is necessary 
to perform a full renormalisation group analysis. Fortunately this is fairly simple since 
the renormalisations of h and A do not mix. This is because the renormalisation of h 
may be discussed in terms of its contribution to the propagator ((f)(x, t)(f)(x', t')), to which 
A does not contribute, while, since A is a bulk coupling, its renormalisation cannot be 
affected by the localised interaction h(p(f)(0). The bulk renormalisation of A is discussed 
in Ref. 3. Here we use the notation l R to denote the dimensionless renormalised coupling 
\RK~ 2e '/ dc ( k \ where e' = d c (k) — d. The renormalisation of h is then carried out for 
A = 0. Since this is a Gaussian theory, this is simple, but non-trivial due to the localised 
form of the interaction. The renormalised coupling h R may be defined in terms of the 
truncated Fourier-Laplace transform of (<j>(x, £)0(O, t f )<j)(0, t')(f)(x", t")), evaluated at the 
normalisation point imaginary frequency s = k 2 . We thus find 

h R = h(l + hI d (K))- 1 

where 

, . f 1 d d p n~ e 

Jd « = / -—^T^d 2 = B ( d ) 

J s + p z (27r)" s=k 2 e 

The dimensionless coupling g R = h^K^ e (where now e = 2 — d) then has the beta function 

Pg(9R) = -^R + B(d)g 2 R 

where 5(2) = l/2n. This is exact to all orders in g R . In addition to the coupling constant 
renormalisation, however, it is necessary to perform a multiplicative renormalisation of 
0(0, t). This may be seen by considering the correlation function (0(O,£)0(p = 0,t = 0)), 
whose Laplace transform is s _1 (l + /i/^s 1 / 2 )) -1 , of which the divergence for d = 2 cannot 
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be removed by the renormalisation of h. We therefore define </>r(0, t) = Zq4>(0, t) to remove 
this factor, where Zq = 1 + (B(d)/e)hK~ e . 

Consider now (t, g R , £ R , k) = (4>n(0,t)). This satisfies a renormalisation group 



equation 



° - T0(9R)) C$ = 



'dt R nwwj^R 
where 70 = (nd/dx) lnZg = —B(d)g R . The solution as t — > 00, for d < 2, is 

C$\t,g R ,i R ,K) ~ (k 2 *)-^ 2 ^)- 6 /^^- 2 , (3) 

where (7* = e/B{d), and is the running coupling The first prefactor on the right 
hand side comes from the canonical scaling dimension of <p, the second from the anomalous 
dimension 70 (<7*)- 

For k > 2, there is a regime where d c (k) < d < 2. In this case, A is irrelevant and 
A ~ t e / d c(k)^ w ith e ' < 0. The correlation function of the right hand side of (3) is then 
given by the sum of tree diagrams, equivalent to solving the inhomogeneous rate equation 
(1) (with (f> replacing n) with the boundary condition K e g*(/)(0) = lim r ^o Sd rd ~^d(f)/dr, 
which comes from varying with respect to 0(0) and integrating by parts. is the area 
of a unit (i-dimensional sphere, and A = £n e . The solution is proportional to A~ dc ( fc )/ 2 , so 
that C R \t) ~ l/t 1+e '/ 2 . As before, the higher cumulants are irrelevant for d > d c (k), so 
we obtain the stretched exponential result for P(t) given in the summary below. At d = 2, 
the prefactor (/t 2 t)~ e / 2 is replaced by (ln(« 2 t)) _1 . This results in an extra (hit) -1 factor 
in the exponent. 

When d < d c (k), I also flows towards a non-trivial fixed point £* = O(e'), so that the 
correlation function on the right hand side of (3) is asymptotically independent of t. The 
prefactors combine to give a simple 1/t dependence, which integrates to give Int. The 
amplitude of this term is given, to leading order in e' by setting £ = I* in the solution of 
(1), which gives a universal amplitude of 0(e' rfc ( fc )/ 2 ). Corrections to this come from loop 
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corrections to the right hand side, and, more importantly, from the higher order cumulants, 

which satisfy similar renormalisation group equations and whose integrals all scale like 

/ iik— 2)/(k— 1) 

Int. However, their amplitudes are suppressed for small e by powers of e y " K ' . For 

d = d c (k), £(t) flows to zero like (In*) -1 , so that ~ (l/t)(lnt) 1 /(fc-l), with the higher 

order cumulants being suppressed only by powers of In t. 

To summarise the different cases when k > 2, we have 
( 3-const.i 1 " MM d> 



P(t) = { 



g-const.i 1 "^^)/ 2 / 111 * d = 2- 
e-™*^ Mk))/ \ 2>d>d c (k); 



e -const.(lnt) fe /( fe ^ d = d c (k); 
> t- 9 , d < d c (k), with = 0(e /_4(fc)/2 ). 

Turning now to = 2, the case d > 2 has already been discussed. For d < 2 the 
same arguments as for k > 2 show that P(t) ~ t - ^, with a universal exponent. However, 
the dependence on e is now of a different form. To leading order, we may solve the rate 
equation to calculate the right hand side of (3), setting 1 = 1* and g* = 0. (Higher order 
terms in g* are higher order in e.) This means that the amplitude is, to leading order, 
that of the bulk density 3 (<j>(t)) ~ l/(4ne)t. In addition, there is a factor of hZ^ 1 ~ 2ne 
in relating h(<j)(0,t)) to C^\t). Thus the factors of e cancel, and we find 

9= 1 - + 0(e) (4) 

The case d = 2 is the most interesting, since both h and A are marginally irrelevant 
there. In this case the prefactor behaves as (47r//i)(ln£) -1 , while 0(0, t) (to leading order 
in gji) behaves as 3 (l/87r)(lnt/t). Thus we get a competition between the two running 
couplings which results in a power behaviour for P(t), with 9 = |, consistent with its 
limit as e — > 0. Although (4) has the appearance of a conventional mean field result with 
0(e) corrections below d = d c = 2, this is not the case: the value of 9 for <i = 2 comes 
about by a subtle cancellation of fluctuation effects, and the exponent, for d > d c , is not 
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universal. In addition, for d = 2, the corrections to C^\t) are suppressed by powers of hit 
only. Thus the leading corrections to lnP(t) are proportional to J"*(l/t' lnt')dt' ~ In Int. 
This will give a logarithmic prefactor multiplying the power law t~ 1//2 . Unfortunately, the 
calculation of the exponent of this logarithm requires a two-loop calculation which is more 
difficult. 

In principle, it is possible to compute higher order terms in the e-expansion (4) . How- 
ever, a similar expansion 3 for the bulk density amplitude does not appear to extrapolate 
well to d = 1. The O(e) corrections to (4), which come from the second cumulant in (2), are 
expected to be negative, consistent with the result of Derrida et a/. 1 , who find that 9 as 0.37 
for d = 1. It is also straightforward to extend the analysis of the case k = 2 to include the 
reaction A + A — > A. If the rate for this second process is A', the effect is to change the 
interaction part of the action S to (2 A + X')^ 2 + (A + \')(p(p 2 . This may be brought back 
to the standard form 3 by rescaling <f> = = where £ = 2(A + A')/(2A + A'). 

The result is that the density amplitude for d < 2, and therefore the exponent 6 to lowest 
order, is modified by this factor of £. 3 In one dimension, the reaction A + A — > A occurs 
in the domain wall dynamics of the (/-state Potts model with q ^ 2, with A'/A = q — 2. 
Thus, to leading order in e, the power is modified to 

e = ^—± + o(e) 

so that, at least to this order, 9 increases with q as found for d = 1. The (/-dependence is, 
however, more complicated for the higher order terms. 
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